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, The large-distance behaviors of the random field and random anisotropy O(N) models are studied 

■ with the functional renormalization group in 4 — e dimensions. The random anisotropy Heisenberg 

(TV = 3) model is found to have a phase with an infinite correlation length at low temperatures and 
weak disorder. The correlation function of the magnetization obeys a power law (m(ri)m(r2)} ~ 
O ' n — r2|~ 0,62E . The magnetic susceptibility diverges at low fields as \ ~ H~ 1+0 15e . In the random 

field O(N) model the correlation length is found to be finite at the arbitrarily weak disorder for any 
JV > 3. The random field case is studied with a new simple method, based on a rigorous inequality. 
This approach allows one to avoid the integration of the functional renormalization group equations. 
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I. INTRODUCTION 



The effect of impurities on the order in condensed matter is interesting, since the disorder is almost inevitably 
present in any system. If the disorder is weak the short-range order is the same as in the pure system. However, 
the large-distance behavior can be strongly modified by the arbitrarily weak disorder. This happens in the systems 
of continuous symmetry in the presence of the random symmetry breaking field |Q. The first experimental example 
of this kind is the amorphous magnet Oil. During the last decade a lot of other related objects were found. These 
are liquid crystals in the porous media H, nematic elastomers ||, He-3 in aerogel || and vortex phases of impure 
i superconductors [Q. The nature of the low-temperature phases of these systems is still unclear. The only reliable 

statement is that a long-range order is absent |Hp|-|lO[]. However, other details of the large- distance behavior are 
C*~) . poorly understood. 

The neutron scattering [pd[| reveals sharp Bragg peaks in impure superconductors at low temperatures and weak 
external magnetic fields. Since the vortices can not form a regular lattice M it is tempting to assume that there 
is a quasi-long-range order (QLRO), that is the correlation length is infinite and correlation functions depend on 
the distance slow. Recent theoretical |l2|,[l3| and numerical JlJ] studies of the random field XY model, which is 
the simplest model of the vortex system in the impure superconductor [Q, support this picture. The theoretical 
Q\ | advances 12,13|] are afforded by two new technical approaches: the functional renormalization group |l5| ] and the 



0^ ■ replica variational method Jl6|. These methods are free from drawbacks of the standard renormalization group and 
give reasonable results. The variational method regards the possibility of spontaneous replica symmetry breaking 
and treats the fluctuations approximately. On the other hand the functional renormalization group provides a subtle 
analysis of the fluctuations about the replica symmetrical ground state. Surprisingly, the methods give close and 
sometimes even the same results. 

Both techniques were originally suggested for the random manifolds pR 16[ and then allowed to obtain information 



about some other disordered systems with the Abelian symmetry [|12 1 3] 1 7 -fi 9[] . Less is known about the non-Abelian 
systems. The simplest of them are the random field (RF) 0] and random anisotropy (RA) Q Heisenberg models. 
The latter was introduced as a model of the amorphous magnet ^,|) . In spite of a long discussion, initiated by Ref. 
[20 1, the question of QLRO in these models is still open. There is an experimental evidence in favor of no QLRO 
[21|. On the other hand recent numerical simulations [^2| support the possibility of QLRO in these systems. The 
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only theoretical approach, developed up to now, is based on the spherical approximation ||^,^4| and predicts the 
absence of QLRO at W > 1 magnetization components. However, there is no reason for this approximation to be 
valid at N ~ 1. 

In this paper we study the RF and RA O(N) models in 4 — e dimensions with the functional renormalization group. 
The large-distance behaviors of the systems are found to be quite different. Whereas in the RF O(N) model with 
N > 3 the correlation length is always finite, the RA Heisenberg (N = 3) model has a phase with QLRO. In that 
phase the correlation function of the magnetization obeys a power law and the magnetic susceptibility diverges at low 
fields. 

The paper has the following structure. In the second section the models are formulated and a qualitative discussion 
is given. The third section contains a derivation of the one-loop renormalization group (RG) equations. The 4th 
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section is devoted to the RF model. The absence of QLRO in that model at N > 3 is shown with a new method, 
based on a rigorous inequality. This approach simplifies tedious RG calculations and can be useful in other problems. 
The RA case is considered in the 5th section. The stable RG fixed point corresponds to QLRO in 4 — e dimensions 
at N ~ 1. In particular, at the weak disorder the correlation length is infinite in the low temperature phases of 
the RA XY (N = 2) and Heisenberg (N — 3) models. However, QLRO is absent at N > 10. In 4 dimensions the 
correlation functions of the RA Heisenberg model depend on the distance logarithmically. The exact result for the 
two-spin correlator is given by the In -0 ' 62 R law. The Conclusion contains a discussion of the results. Appendix A is 



devoted to a generalization of the Schwartz-Soffer inequality 25 1 . The generalized inequality is applied to the stability 



analyses of the RG fixed points. Appendix B describes a simple Migdal-Kadanoff renormalization group approach 
that reproduces qualitatively the results of the rigorous method. This approximation provides good estimations of the 
critical exponents of the RA XY and Heisenberg models. Appendix C includes some technical details of the functional 
RG in the spherical model. 



II. MODEL 



To describe the large-distance behavior at low temperatures we use the classical nonlinear cr-model with the Hamil- 
tonian 

H= /'d D 2 ;[J^9 M n(x)9 M n(x) + Vi mp (x)], (2.1) 

J |U 

where n(x) is the unit vector of the magnetization, T4 mp (x) the random potential. In the RF case it has the form 

^ mp = - M*K(x); a = 1, N, (2.2) 

a 

where the random field h(x) has a Gaussian distribution and (/i ct (x)/i l a(x / )) = A 2 5(s. — x!)S a p. In the RA case the 
random potential is given by the equation 

Mmp = - y j T a p (x) n a (x) ng (x) ; a, ft = 1, ...,N, (2.3) 

where r a ^(x) is a Gaussian random variable, (r Q/ g(x)r 7 ,5(x')) = A 2 5 a j6p$6(yL — x'). The random potential fl2.3| ) 
corresponds to the same symmetry as the more conventional choice V\ mp = — (hn) 2 but is more convenient for the 
further discussion. 

We assume that the temperature is low and the thermal fluctuations are negligible. The Imry-Ma argument [QJ^] 
suggests that in our problem the long-range order is absent at any dimension D < 4. One can estimate the Larkin 
length, up to which there are strong ferromagnetic correlations, with the following qualitative RG approach. Let one 
remove the fast modes and rewrite the Hamiltonian in terms of the block spins, corresponding to the scale L = ba, 
where a is the ultraviolet cut-off, b > 1. Then let one make rescaling such that the Hamiltonian would restore its 
initial form with new constants A(L), J(L). Dimensional analysis provides estimations 

J{L) ~ b D - 2 J{a);A(L) ~ b D/2 A(a) (2.4) 

To estimate the typical angle (j) between neighbor block spins, one notes that the effective field, acting on each spin, 
has two contributions: the exchange contribution and the random one. The exchange contribution of order J(L) is 
oriented along the local average direction of the magnetization. The random contribution of order A(L) may have any 
direction. This allows one to write at low temperatures that <j>(L) ~ A(L)/J(L). The Larkin length corresponds to 
the condition 4>(L) ~ 1 and equals L ~ (J/A) 2 /( 4 ~ D ' in agreement with the Imry-Ma argument [p. If Eq. ( p~i| ) were 
exact th e La rkin length could be interpreted as the correlation length. However, there are two sources of corrections 



to Eq. (2.4). Both of them arc relevant already at the derivation of the RG equation for the pure system in 2 + e 
dimensions pq |. The first source is the renormalization due to the interaction and the second one results from the 
magnetization rescaling which is necessary to ensure the fixed length condition n 2 = 1. The leading corrections to 



Eq. (2.4) are proportional to 4> 2 J,4> 2 A. Thus, the RG equation for the combination (A(L) j J (L)) 2 reads 

4 



d (A{L)\ = (A{L)\ , (ML)Y 
d\nL\j{L)) e \j{L)) +C \J(L)J 



D (2.5) 
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If the constant c in Eq. (2.5) is positive the Larkin length is the correlation length indeed. But if c < the RG 
equation has a fixed point, corresponding to the phase with an infinite correlation length. As seen below, both 
situations are possible, depending on the system. 

The large-distance behaviors of the RF and RA O(N) models are known in two limit cases: N — 2 and N = oo. 
In the spherical limit (N = oo) QLRO is absent (Appendix C, pj|) while the XY model possesses QLRO Jl^ , |f3| , p7t . 
Hence, the ordering disappears at some critical number N c of the magnetization components. This critical number 
is larger in the RA model, since the fluctuations of the magnetization are stronger in the RF case. Indeed, in the RF 
model the magnetization tends to be oriented along the random field, whereas in the RA case there are two preferable 
local magnetization directions so that the spins tend to lie in the same semispace. 



III. RG EQUATIONS 



In the previous section the RG equations are discussed from the qualitative point of view. Eq. (2.5) corresponds 
to the Migdal-Kadanoff approach of Appendix B. In the present section we derive the RG equations in a systematic 
way. 

The one-loop RG equations for the A-component RF and RA models in 4 + e dimensions were already derived 
in Ref. p8| ]. We can directly use that result, since the RG equations in dimensions D < 4 can be obtained by just 
changing the sign of e. However, the approac h P8fl is cumbersome and we provide below a simpler derivation. 

We use the method, suggested by Polyakov pa] for the pure system in 2 + e dimensions. This method is technically 
simpler and closer to the intuition than the other approaches. A disadvantage of the method is the difficulty of the 
generalization for the higher orders in e. This generalization requires the field-theoretical approach p9| . 

The same consideration as in the XY jl3| and random manifold [[l5) models suggests that near a zero-temperature 
fixed point in 4 — e dimensions there is an infinite set of relevant operators. Let us show that after the replica averaging 
the relevant part of the effective replica Hamiltonian can be represented in the form 

h r = f d D x[j2 ^ E - E %A (3-1) 

a fl ab 

where a, b are replica indices, R{z) is some function, T the temperature. We ascribe to the field n the scaling 



dimension 0. We also assume that the coefficient before the gradient term in (3.1) is 1/(2T) at any scale. Then in the 
(4 — e)-dimensional space the scaling dimension of the temperature = — 2 + 0(e). Any operator A containing m 
different replica indices is proportional pl| to 1/T m . Hence, the scaling dimension of the operator A satisfies the 
relation A^ = 4 — n + mAj + 0(e), where n is the number of the spatial derivatives in the operator. The relevant 
operators have A^ > 0. Hence, the relevant operators can not contain more than two diff eren t replica indices. A 
symmetry consideration shows that all the possible relevant operators are included into Eq. ( |3.l| ). The function R(z) 
is arbitrary in the RF case. In the RA case R(z) is even due to the symmetry with respect to changing the sign of 
the magnetization. 

The one-loop RG equations for the A-componcnt model in 4 — e dimensions are obtained by a straightforward 
combination of the methods of Refs. |l5j and |2(|. We express each replica n a (x) of the magnetization as a combination 
of fast fields </>"(x), i — 1, N — 1 and a slow field n' a (x) of the unit length. We use the representation 




n»(x) = n'«(x) l-^(x)) 2 l^^(x)e«(x), (3.2) 



where the unit vectors e"(x) are perpendicular to each other and the vector n /a (x). The slow field n /a can be chosen 
in different ways. For example, one can cut the system into blocks of size L a, where a is the ultra-violet cut-off. 
In the center x of a block the vector n' a (x) should be parallel to the total magnetization of the block. In the other 
points the field n' a should be interpolated. We assume that the fluctuations of the magnetization are weak, that is 
(</>?) <C 1. Then the fluctuations of the field n a are orthogonal to the vector n' a because of the fixed length constraint 
(n a ) 2 = 1. 



To substitute the expression (3.2) into the Hamiltonian we have to differentiate the vectors ef. We use the following 
definition 
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It is easy to show |26|] that ^^(c^) = X^(^M n ' a ) ■ With the accuracy up to the second order in <p the replica 



Hamiltonian (3.1) can be represented as follows 

h r = f d D x[^ £{(^n'°) 2 (i (^) 2 ) + + - ftMl?) 

~ ]T{i?(n'V b ) + A a \cj>tf + Btftitf + C$ </>?$}], (3.4) 

ab 

where the summation over the repeated indices i, j, k, \i is assumed and 

A ab = _( n 'V 6 )i2V° n ' 6 ); S |f = (n' b <)(n' b e;*)i?''(n' a n' b ); 

C?J> = (e^)R'(n' a n' b ) + (n' Q e b J(n' b <)i?''(n' a n' b ). (3.5) 

In the last formula R' and R" deno te the first and second derivatives of the function R(z). We have omitted the terms 
of the first order in <f> in Eq. ( |3.4| ) . These terms are proportional to the products of the fast field (f> and some slow 
fields. Hence, they are non-zero only in narrow shells of the momentum space. One can show that their contributions 
to the RG equations are negligible. 

To obtain the RG equations we have to integrate out the fast variables </>" . Near a zero-temperature fixed point the 
Jacobian of the transformation n — ► (n', </>,) can be ignored, since the Jacobian does not depend on the temperature. 
The integration measure is determined from the condition that t he f ield n /a is a slow part of the magnetization. Th is 



condition imposes restrictions on the fields <j>. The expression ( |3.4| ) depends on the choice of the vectors ef (p.2j). 
However, after integrating out the fields <fi the Hamiltonian can depend only on the slow part n /a of the magnetization. 
One can make the calculations simpler, considering special realizations of the field n'°. To find the renormalization 



of the disorder-induced term R(z) (3.1) we can assume that the field n' a does not depend on the coordinates. The 



renormalization of the gradient energy can be determined, assuming that the vectors n' a (x) depend on o ne spatial 



coordinate only and lie in the same plane. In both cases the vectors ef can be chosen such that in Eq. ( p.3[) /" ik = 0. 
At such a choice the integration measure can be omitted and the fields 4>f can be considered as weakly interacting 
fields with the wave vectors from the interval 1/a > q > 1/L. 

To derive the one-loop RG equations we express the free energy via the Hamiltonian ( |3.4| ). Then we expand 
the exponent in the partition function up to the second order in (Hr — f d D xJ2^ i (d fJ ,4>i) 2 /(2T)) and integrate over 
cj>f. Finally, we make a rescaling. The vectors ef can be excluded from the final expressions with the relation 
Sj(ae")(be") = (ab) — (an' a )(bn' a ). In a zero-temperature fixed point the RG equations are 

-^j- = -(D - 2) + 2(N - 2)R'(1) + 0(R 2 ,T); (3.6) 

= = eR(z) + 4(N - 2)R(z)R'(l) - 2(N - l)zR'(l)R'(z) + 2(1 - z 2 )R' (1)R" (z) 

cl In 1j 



+ (R'(z)f(N - 2 + z 1 ) - 2R'(z)R"(z)z(l - z 2 ) + (R"(z)) 2 (l - z 2 ) 2 , (3.7) 



where the factor l/(87r 2 ) is absorbed into R(z) to simplify notations. The RG equations become simpler after the 
substitution of the argument of the function R(z): z = cos<f>. In terms of this new variable one has to find even 
periodic solutions R((j>). The period is 2ir in the RF case and tt in the RA case due to the additional symmetry of the 
RA model. The one-loop equations get the form 

d In T 

-j^— = -(D - 2) - 2(N - 2)R"(0) + 0(R 2 ,T); (3.8) 
d In 1j 

= = eR(<fi) + (R"(<i>)) 2 - 2R"(<t>)R"(Q) - 

(N - 2)[4R((j>)R" (Q) + 2ctg<t>R'(<t>)R"(0) - } + 0(R 3 , T) (3.9) 



Eq. (3.8) provides the following result for the scaling dimension At of the temperature 

A T = -2 + e-2(A r -2) J R"(0). (3.10) 
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The two-spin correlation function is given in the one-loop order pa] by the expression 



(n a (x)n a (x')> = (n' a (x)n' a (x'))(l 



(#) 2 »- 



(3.11) 



Hence, in the fixed point (n(x)n(x')) ~ |x — x'| T> , where 

r} = -2(N - l)R"((p = 0) 



(3.12) 



Let us hnd the magnetic susceptibility in the weak uniform external field H. We add to the Hamiltonian (3.1) the 
term — f d D xHn^/T (the field is dir ecte d along the z-axis). The renormalization of the field H is determined by 
the renormalization of the temperature (3.S) and the field n. In the zero-loop order the renormalized magnetic field 
h(L) depends on the scale as h(L) = H X (L/a) 2 . Hence, the correlation length R c ~ H~ x / 2 . The magnetization, 
averaged over a block of size i? c , is oriented along the field. The absolute value of this average magnetization is 
proportional to Rc V ^ 2 ■ This allows us to calculate the critical exponent 7 of the magnetic susceptibility x(H) ~ H ~ 7 
in a zero-temperature fixed point: 



7 = 1 + (N - 1)R"{4> = 0)/2. 



(3.13) 



In Ref. [M Eqs. (3. £,3. 7) were derived with a different method. In that paper the critical behavior in 4 + e 
dimensions was studied by considering analytical fixed point solutions R(z). In the Heisenberg model, analytical 
solutions are absent and they are unphysical for N 3 p8[ . In 4 — e dimensions appropriate analytical solutions are 
absent for any N. To demonstrate this let us differentiate Eq. (3.7) over z at z — 1. For any analytical R{z) we 
obtain the following flow equation 



dR'(z = 1) 
d\nL 



eR'(z= 1) + 2{N -2)(R'{z = l jf 



(3.14) 



At N > 2 th e fix ed point of this equation R'(z = 1) = — e/[2(N — 2)] < 0. It corresponds to the negative critical 
exponent 77 ( |3.12 ) and hence is unphysical. However, we shall see that in the RA model some appropriate non- 
analytical fixed points R(z) appear. In these fixed points R"(z — 1) = 00. In Ref. p8| the RG charges are the 
derivatives of the function R(z) at z — 1. Thus, in a non-analytical fixed point these charges diverge. In the systems 
with a finite number of the charges their divergence implies the absence of a fixed point. However, if the number of 
the RG charges is infinite such a criterion does not work and is even ambiguous. Indeed, the set of charges can be 
chosen in different ways and e.g. the coefficients of the Taylor expansion about z = remain finite in our problem. 



IV. RANDOM FIELD 

For the RF XY model the one- loop RG equations ( |3.8| , |3.9| ) can be solved exactly The solution corresponds 

to QLRO with the critical exponents r\ = ir 2 /9e,j = 1 — 7r 2 /18e. In the first order in e the exponent r\ equals the 
prefactor C before the logarithm in the correlation function |13| of the angles ^(x) between the spins n(x) and some 
fixed direction: ((</>(xi) — </>(x2)) 2 ) = Cln|xi — xa|. We expect that this coincidence does not extend to the higher 
orders. 



If N 7^ 2 the RG equation (3.9) is more complex. Fortunately, at N > 3 there is still a simple method to study the 
large-distance behavior. The method is based on the Schwartz-Soffer inequality [g5| and shows that QLRO is absent. 

In Ref. [p5| the inequality is proven for the Gaussian distribution of the random field. It can also be proved for the 
arbitrary RF distribution (Appendix A). 

Let us demonstrate the absence of physically acceptable fixed points in the RF case at N > 3. We derive some 
inequality for critical exponents. Then we show that the inequality has no solutions. We use a rigorous inequality for 
the connected and disconnected correlation functions |2q] 

((n(q)n(-q))) = (n (q)n (-q)) - (n (q)n 6 (-q)) < const yj (n a (q)n a (-q)), (4.1) 

where n(q) is a Fourier-component of t he m agnetization, a, b are replica indices. The disconnected correlation function 
is described by the critical exponent ( 3.1 2j ) . The large- distance behavior of the connected correlation function in a 



zero-temperature fixed point can be derived from the expression 

X ~ [{(n(0)n(x)))d D x (4.2) 
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and the critical exponent of the susceptibility ( 3.1 3| ) . The integral in the right hand side of Eq. (4.2) is proportional to 
R®~ Vl , where R c is the correlation length in the e xtern al field H , 7/1 the critical exponent of the connected correlation 
function. For the calculation of the exponent 7 ( 3.1 3| ) we used the zero-loop expression of R c via H. Now we need 
the one-loop accuracy. In this order R c ~ iJ^ 1 /[ 2 ~( Ar ~ 3 )- R (°)]. This allows us to get the following equation for the 
exponent 771 



D — 2 — 2R"(0). 



(4.3) 



In a fixed point Eq. (4.1) provides an inequality for the critical exponents of the connected and disconnected correlation 
functions p5|. The inequality has the form 



2(2- D + m) >4-D + V . 



This allow us to obtain the following relation 



3- N 

i-D<j^r, + o(R), 



(4.4) 



(4.5) 



where r\ is given by Eq. ( 3.12j ). The two-spin correlation function can not increase up to the i nfin ity as the distance 
increases. Hence, the critical exponent r\ is positive. At N > 3 this is incompatible with Eq. (4.5) at small e. Thus, 
there are no accessible fixed points for N > 3. This suggests the strong coupling regime with a presumably finite 
correlation length. 



Certainly, in the RF XY model |12| , |l3f Eq. (4.5) is satisfied. However, the unstable fixed points of the RG equations 
fL3| do not satisfy the inequality. 

The marginal Heisenberg case N = 3 is the most difficult, since in the one-loop order the right hand side of Eq. 
( [4.5P equals zero at N = 3. Hence, the two-loop corrections may be relevant. The RF Heisenberg model is beyond 
the scope of the present paper. 



V. RANDOM ANISOTROPY 



In this section we investigate the possibility of QLRO in the RA O(N) model. The first subsection is devoted to 
the simplest case of the XY model. The second sub section contains an inequality for the critical exponent r/. The 
derivation of the inequality is analogous to Eq. (4.5). This inequality is applied in the next subsections. The third 
subsection contains the results for the Heisenberg model. In the last subsection we consider the case N > 3. 



A. N 



This case is studied analogously to the RF XY model At N = 2 the RG equation (3.S) can be solved 

analytically. Its solution is a periodical function with period it. In interval < <p < ir the fixed point solution R(<f>) is 
given by the formula 



Ri4>) = TU V36- (^(I-^/tt)) 1 



(5.1) 



It is a stable fixed point. This can be verified with the linearization of the flow equation fl3.9| ) for the small deviations 
from the fixed point. Another proof of the stability is based on the inequality of the next subsection. 

The stable fixed point corresponds to the QLRO phase at low temperatures and weak disorder. The critical 
exponents i] = 7r 2 e/36, 7=1 — ir 2 e/72. 

The solution (5.1) is non-analytical at <fi = 0, since R IV ((/> = 0) = 00. Hence, the Taylor expansion over <f> is absent. 
However, a power expansion over \(j)\ exists. We shall see below that the same behavior at small <j> conserves also at 
other N. 



B. An inequality for a critical exponent 



We use the same approach as in the RF model. Since in the RA case the random field is conjugated with a second 
order polynomial of the magnetization, the Schwartz-Soffer inequality pa] should be applied to correlation functions 



G 



of the field m(x) = (n 2 (x)) 2 — 1/N, where n z denotes one of the magnetization components, 1/N is subtracted to 
ensure the relation (to) = 0. 

To calculate the critical exponent \x of the disconnected correlation function we use the representation (3.2) and 
obtain the relation 

<m (x)m<V)> = (m'°(x)m"V)) (l - ^^f^ ) > ( 5 - 2 ) 

where a is a replica index, to' — (n' z ) 2 — 1/N the slow part of the field to. One finds /i = —4NR"(0). 

The critical exponent /ii of the connected correlation function is determined analogously to the RF case. We apply 
a weak uniform field H , conjugated with the field to, and calculate the susceptibility dm/ dH in two ways. The result 
for the critical exponent is \i x = D - 2 - 2(N + 2)R"(0). 

Th e S chwartz-Soffer inequality provides a relation between the exponents /i and fi±. It has the same structure as 
Eq. (4.4). Finally, we obtain the following equation 

V>^-(N-l)+o(R). (5.3) 

In terms of the RG charge R(4>) this inequality can be rewritten in the form 

R"(0) < -e/8 + o(R). (5.4) 



C. N = 3 



In this case we solve Eq. ( |3.9| ) numerically Since coefficients of Eq. (^9) are large as <\> — > 0, it is convenient to use 
a series expansion of the fixed-point solution R(4>) at small <fi. At the larger <p the equation can be integrated with 
the Runge-Kutta method. The following expansion over t = (1 — z)/2 = | sin(^/2)| holds 



(N — l)a 2 
l-4(JV-2)o 



2a sin — ± 
2 

2a 

y 



4^2 l-a + 2(N -2)a 2 



3(7V + 4) 



N ■ 



•3^1 
Sin — 
2 1 



0(|sin 5 ||), 



(5.5) 



Similar to the random manifold 



where a — R"(<fi = 0)/e. We see that the RG charge R(4>) is non-analytical at small < 
H and random field XY ||] models R IV (0) = oo. 

Numerical calculations sho w th at at any N the solutions, compatible with the inequality (5.4), have sign "+" 
before the third term of Eq. ( |5.5| ). The solutions to be found are even periodical functions w ith period n. Hence, 
their derivative is zero at0 = 7r/2. AtiV = 3 there is onl y one solu tion that satisfies Eq. (|]J). It corresponds 
to R"{4> = 0) = — 0.1543e. If this solution is stable Eqs. ( 3.12| ,3.13) provide the following results for the critical 
exponents 



T) = 0.62e;7 = 1 - 0.15e. 



(5.6) 



All the other solutions of Eq. (3.9) do not satisfy Eq. (5.4) and hence are unstable. 

We have still to test the stability of the solution found. For this aim we use an approximate method. First, we find 
an approximate analytical solution of Eq. ( |3.9| ). We rewrite Eq. ( |3.9D , substituting uj(R"((/))) 2 for (R"((f>)) 2 . The case 
of interest is u> = 1 but at u> = the equation can be solved exactly. The solution at u) = 1 can then be found with 
the perturbation theory over u>. The exact solution at u> = is Ruj=o(4>) — e(cos 2<^>/24 + 1/120). The corrections of 
order oj k are trigonometric polynomials of order 2(k + 1). The first correction is 



_ , 2uj€ toe 
R\{<p) = — — — - cos 2(f>+ —— cos 40 -I- const 
99 264 



(5.7) 



After the calculation of the corrections we can write an asymptotic series for the critical exponent rj ( 3.1 2[ ) : rj = 
e(0.67 — 0.08w + O.lAu; 2 — ...). The resulting estimation ij = e(0.67 ± 0.08) agrees with the numerical result (5.6) 
well. This allows us to expect that the stability anal ysis of the solution i? w= o of the equation with u> = provides 
information about the stability of the solution of Eq. (3.9). 
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To study the stability of the exact solution of the equation with u) = is a simple problem. We introduce a small 
deviation r (</>): R(4>) = Ru}=o{<P) + f(4>) and write the flow equation for this deviation: 



rflnL 



= (br{4>) +r"(4>) + r"(0) cos 2</>)/3 + const x r"(0). 



(5.8) 



It is convenient to use the Fourier expansion r(<ft) = ^2 m a m cos2m(f). The flow equations for the Fourier harmonics 
can be easily integrated. We see that a m — > as L — > oo for any m > 0. The solution is unstable with respect to 
the constant shift ciq. However, this instability has no interest for us, since the correlation functions do not change 
at such shifts |l5). Indeed, the constant shift corresponds to the addition of just a random term, independent of 
the magnetization, to the Hamiltonian (2.1). Thus, the RG equation possesses a stable fixed point. This fixed point 
describes the QLRO phase with the critical exponents (5.6). 

In the Abelian systems the results of the functional RG are supported by the variational method p6| . In our 
problem this method can not be applied. However, it is interesting that in the Abelian systems the functional RG 
equations without (R"(<fi)) 2 reproduce the variational results. 

As usual in critical phenomena, in 4 dimensions the one-loop RG equations allow one to obtain the exact large- 
distance asymptotics of the correlation function. In the 4-dimensional case R(<f) = R(4>) / In L, where R(<fi) satisfies 
Eq. (3.9) at e = 1. We obtain the following result for the two-spin correlation function with Eq. (3.11) 



(n(x)n(x')) 



ln-°- 62 lx 



(5.9) 



D. N > 3 

Numerical analysis of Eq. (|3.9|) shows that solutions, compatible with Eq. (5.4), are absent at N > 10. Hence, 



QLRO is absent for any N > 10. In the spherical model (N — oo) the absence of fixed points can be demonstrated 
analytically (Appendix C). This agrees with the previous results P^3j - For each integer N < 10 the RG equation 
( p.9[ ) has exactly one solution, satisfying the inequality of section |V B[ These solutions are described in Table I. In 
t he ta ble, rj is the critical exponent of the two-spin correlation function, At the scaling dimension of the temperature 

( B) . 



Unfortunately, it is not clear if the fixed points, found at TV > 3, survive in 3 dimensions. A zero-temperature fixed 
point can exist only if the scaling dimension of the temperature is negative. Table | shows that scaling dimension is 
positive in the one-loop approximation at e = 1 and N > 5. In the 3-dimcnsional 0(4) model the one-loop correction 
to the scaling dimension — 2(iV — 2)i?"(0) 0.7e is close to the zero-loop approximation —2 + e. Thus, the next orders 
of the perturbation theory are crucial to understand what happens in 3 dimensions. 



In the 0(2) model the scaling dimension At = — 2 + e is exact J_3 15 1. Hence, QLRO disappears in 2 dimensions. 
In systems with a larger numbers of magnetization components fluctuations become stronger. Thus, one expects the 
absence of QLRO in all the two-dimensional O(N) models. 

At the zero temperature Eq. ( |3.9| ) is valid independently of the scaling dimension At- It is tempting to assume 
that at the zero temperature QLRO still exists in the RA 0(N > 3) models below the critical dimension, in which 
Ay = 0. However, the experience of the two-dimensional RF XY model does not support such an expectation. Recent 
numerical simulations show that QLRO is absent even in the ground state of that model po| . 

TABLE I. Critical exponents of the RA O(N) model. 



N 


2 


3 


4 


5 


6 7 


8 


9 


V 


ttV36 


0.62e 


Lie 


1.7e 


2.7e 4.6e 


9.0e 


33e 


A T 


-2+e 


-2+1.3e 


-2+1.7e 


-2+2. 3e 


-2+3. 2e -2+4.8e 


-2+8.7e 


-2+30e 
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VI. CONCLUSION 



We have obtained QLRO in the RA Heisenberg model. This is the first example of QLRO in a non-Abelian system. 
The RF disorder tends to destroy the ordering which exists in the RA case. This difference between the RF and RA 
models is not surprising, since the same difference was already obtained in Ref. jl7| for the two-dimensional RF and 
RA XY models with the dipole forces. 

We have not yet dis cussed the role of the topological defects. The contribution of the topological excitations to the 
RG equations (3.8 3~J) is determined by the rare regions where the random field is sufficiently strong to compensate 
the core energy. Hence, similar to the pure system in 2 + e dimensions they are responsible for the non-perturbative 
corrections of order exp(— 1/e). Thus, their effect is negligible at small e. Several studies were devoted to the role of 
the vortices in the RF XY model |3lJ . The theoretical prediction of QLRO in this system is based on the vortexless 
version of the model |l2|,|l3| . A qualitative estimation and variational calculations |H| suggest that the topological 
defects do not change the behavior of the RF XY model at the weak disorder. Our approach allows us to consider 
the XY model, including vortices. We see that QLRO does exist in the model with the defects. 

However, in our problem there may be a more important source of the non-perturbative corrections. The effect 
of the multiple energy minima can lead to corrections of order e 5 / 2 to the RG equations Unfortunately, the 

non-perturbative effects in the RF systems are not well understood. 

The present paper uses a systematic RG approach. However, some results can be reproduced more simply with an 
approximate Migdal-Kadanoff rcnormalization group (Appendix B). 

The question of the large-distance behavior of the RF and RA Heisenberg models was discussed in Ref. [0 on 
the basis of an approximate equation of state. In that paper QLRO was also obtained in the RA case. However, we 
believe that this is an accidental coincidence, since the equation of state |^p| is valid only in the first order in the 
strength of the disorder, while higher orders are crucial for critical properties p4j . In particular, the approach p(J 
incorrectly predicts the absence of QLRO in the RF XY model and its presence in the exactly solvable RA spherical 
model. It also provides incorrect critical exponents in the Heisenberg case. The reason of the mistakes is the fact that 
in the weak external uniform field the perturbation parameter of Ref. pC| ] is large. 

The RA Heisenberg model is relevant for the amorphous magnets 07 At the same time, for their large-distance 
behavior the dipole interaction may be important Besides, a weak nonrandom anisotropy is inevitably present 
due to mechanical stresses. 

In conclusion, we have found that the random anisotropy Heisenberg model has an infinite correlation length and 
a power dependence of the correlation function of the magnetization on the distance at low temperatures and weak 
disorder in 4 — e dimensions. On the other hand, the correlation length of the random field 0(N > 3) model is always 
finite. 
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APPENDIX A: INEQUALITY FOR CORRELATION FUNCTIONS 



In this appendix we derive an inequality for the correlation functions of the disordered systems. We consider the 
system with the Hamiltonian 

H = J da^fffi^x)) - ft(x)m(0(x))], (Al) 

where <fi is the order parameter, h the random field with short range correlations, Hi may depend on some other 
random fields. We prove an inequality for the Fourier components of the field m: 

G co «(q) < const y/Gdisiq), (A2) 



where G<ji s (q) = (m(q)m(— q)) , G con (q) = (m(q)m(— q)) — (m(q))(m(— q)), the angular brackets denote the thermal 
averaging, the bar denotes the disorder averaging. 
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This inequality can be easily obtained in the case of the Gaussian distribution P(h) of the field h p5| . Indeed, in 
the Gaussian case 



Grfis(q) 



P(h) 



dh{q) 



m q (/i) D{h} 



dh(q) ' 



P(/i)m q (/i) D{h} = const / (P(h)h(-q) mci (h)) D{h} 



(A3) 



where J D{h} denotes the integration over the realizations of the random fi eld, m q (/i) = f D{(f>} exp(—H/T)m(q) / 



J D{4>} exp(—H/T). Applying the Cauchy-Bunyakovsky inequality to Eq. (A3) one gets Eq. (A2). 

However, the assumption about the Gaussian distribution of the random field is not necessary. The inequality 
( [A2] ) can also be extended to a more general situation, corresponding to the effective replica Hamiltonian 
Inde ed, i f one adds to any Hamiltonian a weak Gaussian random field h , conjugated with the field m, it suffices for 
Eq. (A2) to become valid. The addition of the Gaussian random field corresponds to the transformation i?(n a n;,) — > 
i?(n a n&) + An a rifc in Eq. (3.1) where A ~ h 2 is a positive constant. Thus, Eq. (A2) is invalid only, if for the arbitrarily 
small A the replica Hamiltonians can not contain the two-replica contribution i?(n a n;,) = i?(n a rifc) — An a ri{,. This 
corresponds to the border of the region of the possible Hamiltonians and has zero probability. 

For systems in the critical domain there is a simple way to understand why the inequality is valid not only in the 
Gaussian case but also in the general situation. This is just a consequence of the universality. 



APPENDIX B: MIGDAL-KADANOFF RENORMALIZATION GROUP 



This appendix contains a simple approximate version of the renormalization group. The results for the critical 
exponents of the XY and Heisenberg models have a very good accuracy. The value of the magnetization component 
number N c , at which QLRO disappears in the RF model, is probably exact. However, the critical number of the 
components in the RA model is underestimated. 



1. Random field 



We use the following ansatz for the disorder-induced term in the Hamiltonian (|2.lD : i?(n Q rifc) = an a nb + 0, where a 
and (3 are constants. This expression corresponds to the Gaussian RF disorder ( |2.2| ) . Below we ignore the generation 
of the other contributions to the function R(z). The mi ssed contributions are related with random anisotropies of 
different orders. In terms of the angle variable 4> (3~§|,3.9) 



R{<j>) = acos0 + (3. 



(Bl) 



To ensure consistency we have to truncate the RG equation (3.9) 
retain only the terms, proportional to cos</> or independent of cf> 
constant a (Bl) 



We substitute the ansatz (Bl) into Eq. (3.9) and 
This leads to the following RG equation for the 



da 
d\nL 



ea + 2a 2 {N -3). 



For N < 3 Eq. (B2) has a stable solution a = e/[2(3 — N)]. The critical exponent (3.12) equals 



V 



(N - l)e 
(3-AO ' 



(B2) 



(B3) 



At N — 2 this result has less than ten percent difference with the systematic theory jl3|. QLRO disappears at N = 3. 
This is the same critical number which is found in section [V. 

For N > 3 a fixed point exists in 4 + e dimensions. It describes the transition between the ferromagnetic and 
paramagnetic phases. In this fixed point the critical exponent (B3) satisfies the modified dimensional reduction 
hypothesis |32|. However, we believe that this is an artifact of the Migdal-Kadanoff approximation, since the correct 
value of the critical exponent differs form Eq. (|B3|). 
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2. Random anisotropy 



In this case we use the ansatz fl(n a n^ = A(n a nt,) 2 + B. In terms of the variable <j) ( pTg|]3.9D R(4>) — a cos 2cf> + (3. 
We again substitute our ansatz into Eq. (3.9) and retain the terms, proportional to cos 20, and the terms, independent 
of <j>. The RG equation for the constant a has the form 



da 
d\nL 



8(iV-6)a 2 



(B4) 



The fixed point solution of this equation is a = e/[8(6 — N)]. It describes the QLRO phase at N < 6. At N — 3 the 
function R((f>) — a cos 2(f> + (3 is just i? w= o of section V C[ The critical exponent of the two-spin correlation function 
is given by the following equation 



e(N-l) 
6-N ' 



(B5) 



At N — 2, 3 this value is close to the results of the systematic approach (Table I). 



APPENDIX C: SPHERICAL MODEL 



In this appendix we consider the spherical RA model with the functional RG. We show that QLRO is absent in 
this model. In the spherical li mit N = oo only the terms, proportional to N, and the term eR(z) should be retained 
in the right hand side of Eq. (3/7). After the change of the variable R(z) = er(z)/N one obtains 

= r{z)[l + 4r'(l)] - 2zr'{l)r'{z) + {r'{z)f. (CI) 

It is convenient to differentiate Eq. ( pl[ ) over z. One gets 

= r'{z)[l + 2r'(l)] + 2r"(z)[r'(z) - zr'(l)}. (C2) 

Analytical functions r(z) can satisfy Eq. (|C2|) at z = 1 only if r'(l) = or r'(l) = —1/2. In both cases Eq. (|C2| ) can 
be easily solved. There are three analytical non-zero solutions: r(z) = — z/2 + 1/4; r(z) = —(1 — z) 2 /4; r(z) — — z 2 /4. 
The last solution only has the necessary symmetry. 

The non-analytical solutions are absent. Indeed, Eq. ( |C2| ) can be integrated with the substitution r'(z) = zt(z). 
The general integral has the form 

W^'"* -a,. (OS) 



(2i(z) + l) 2 *(D+i 

Besides, there are special solutions. They all satisfy the relation t(z) = t(l). Hence, the special solutions are analytical. 
Thus, the function t(z) can be non-analytical at z = 1 only under the condition that z = 1 is a peculiar point of Eq. 
(|C3|). This means that t(l) = or t(l) = —1/2. However, it is easy to verify that in both cases the solution is one of 
the found above. 

We see that the only fixed point of the spherical RA model is R(z) = —ez 2 /(AN). With Eq. (3.12) one finds the 
critical exponent rj = —e/2. Since rj > the solution found is applicable at D > 4. At D < 4 the fixed points are 
absent. Thus, QLRO is absent too. 
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